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Abstract 

For every A-qubit density matrix written in the computational basis, an associated “X-density matrix” can 
be obtained by vanishing all entries out of the main- and anti-diagonals. It is very simple to compute the 
genuine multipartite (GM) concurrence of this associated A-qubit X-state, which, moreover, lower bounds 
the GM-concurrence of the original (non-X) state. In this paper, we rely on these facts to introduce and 
benchmark a heuristic for estimating the GM-concurrence of an arbitrary multiqubit mixed state. By 
explicitly considering two classes of mixed states, we illustrate that our estimates are usually very close 
to the standard lower bound on the GM-concurrence, being significantly easier to compute. In addition, 
while evaluating the performance of our proposed heuristic, we provide the first characterization of GM- 
entanglement in the steady states of the driven Dicke model at zero temperature. 
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1. Introduction 

During the last one and a half decades, increasing interest has been manifested in the topic of multipartite 
entanglement. Just as we learnt, in the early days of quantum information, that bipartite entanglement is an 
invaluable resource for quantum cryptography [T] , communication [5] and speed-up of classical algorithms 131 - 
|B] (to name but a few), we have now a growing awareness of the role that multipartite entanglement plays in 
quantum computing CHin], high-precision metrology mHis], quantum phase transitions |16H25j and even 
biology [MlllZ]- 

Despite significant interest and progresses, many limitations from the theory of bipartite entanglement 
are naturally inherited in the multipartite setting. For example, since the entanglement detection problem 
was shown to be NP-hard already in the bipartite case [2H1 HI] , there is very little hope that computable 
measures of (multipartite) entanglement can be found for an arbitrary quantum state. As a result, much 
of the effort in this field is centered around devising sufficient criteria for (multipartite) entanglement 
detection [301133) . which are then turned into computable (multipartite) entanglement estimators [341137) . 

Amongst all types of multipartite entanglement, special interest is devoted to the detection and quan¬ 
tification of genuine multipartite entanglement, a type of multipartite entanglement established collectively 
between all A parties of an A-partite system. More precisely, a GM-entangled state is any state that cannot 
be written as a mixture of biseparable states, which, in turn, are those states that are separable with respect 
to some bipartition of the relevant Hilbert space [38l |39| . Besides, it is GM-entanglement that occupies the 
highest position in the hierarchy of multipartite entanglement [40) . 
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In this paper, we consider a particular set of sufficient conditions for detecting GM-entanglement intro¬ 
duced by Huber et al. |33j — and turned into an entanglement estimator by Ma et al. |37) — to obtain 
our main result: a heuristic approach for approximating these estimates with a considerably lower compu¬ 
tational cost than that required by the optimization problem appearing in Refs. [31 [37]. In practice, this 
result opens the way for systematical and analytical estimation of GM-entanglement in symmetric multi¬ 
qubit mixed states and, more generally, enables a significant reduction in the computational time required 
by standard numerical solvers, usually accompanied by some small impact on the accuracy of the estimate. 

Our heuristic relies upon three basic facts: (i) the estimates proposed by Ma et al. |37| are actually lower 
bounds on the GM-concurrence — a GM-entanglement monotone introduced by the authors to generalize 
the Wootters concurrence m to the multipartite setting; (ii) this lower bound is saturated for the family 
of fV-qubit X-states, in which case it is also remarkably easy to compute [31] ; and (hi) in general, it 
is computationally easier to determine a local-unitary (LU) transformation that approximates a generic 
density matrix to the X-form, than it is to solve the optimization problem in Refs. [MIET]. Accordingly, 
our so-called X-heuristic, estimates the GM-concurrence of an arbitrary A^-qubit density matrix by applying 
a LU-transformation that minimizes the entries out of its main- and anti-diagonals and, subsequently, 
outputting the GM-concurrence of the associated X-density matrix (obtained by neglecting any entry that 
may have remained out of the main- and anti-diagonals of the LU-transformed density matrix). 

In order to illustrate the strengths and drawbacks of the X-heuristic, we explicitly apply it to the family 
of diagonal symmetric states [431145] and steady states of the driven Dicke model [46H49]. In both cases, we 
discuss the trade-off between computational efficiency and accuracy vis-a-vis the standard scheme for GM- 
concurrence estimation [33]|37]. Our comparative analysis reveals the X-heuristic as a useful computational 
tool for GM-concurrence estimation in multiqubit systems. 

As side results, while assessing the X-heuristic with the steady states of the driven Dicke model at zero 
temperature, we provide the first characterization of multipartite entanglement in this model. In particular, 
we note that our estimates agree on that GM-entanglement is maximal for the parameter values correspond¬ 
ing roughly to a bifurcation of a fixed point in the corresponding semiclassical dynamics, reinforcing the 
long-standing (and widely-accepted) conjecture that multipartite entanglement must be maximal at quan¬ 
tum phase transitions m- While the driven Dicke model is a zero-dimensional many-body problem, the 
bifurcation is controlled by a parameter in the Hamiltonian and the correlations in the quantum steady 
state change at a particular value of that parameter in analogy to the change in the character of a ground 
state in a spatially extended many-body quantum phase transition. Moreover, we show that the two-qubit 
steady states of this model are examples of non-X states whose concurrence can be exactly computed with 
the simple concurrence formula for X-states, settling an open problem posed in Ref m- 

Our paper is structured as follows. In Sec. [^ we review some key concepts relating GM-entanglement 
and its quantification. In Sec. |^our heuristic for GM-concurrence estimation is presented and, in Sec. [^ 
its accuracy is benchmarked against the standard scheme for the two aforemetioned families of A^-qubit 
mixed states. In Sec. [^ we focus on the computational advantage of our heuristic over the standard scheme. 
Finally, in Sec. [^ we summarize our results and discuss some possible avenues for future work. 

2. Preliminaries 

In this section we briefly review the concept of GM-entanglement and some key results regarding its 
detection and quantification. In particular, we outline the scheme introduced in Refs. [SS] |37] (here referred 
to as the “optimal-1$)” method), a current standard for GM-concurrence estimation in multipartite mixed 
states. Although we are only interested in multiqubit states, throughout most of this section we leave the 
dimensionality of each party arbitrary. 

2.1. Genuine Multipartite Entanglement 

An A-partite pure state \ijj) G % = ■ ■ ■®'Hn (with dim.'Hn = d„), is said to be GM-entangled 

if it is not biseparable; \il}) is said to be biseparable, in turn, if there is a bipartition of "H = Ha ®'Hb and 
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a pair of states |^/>^) S 'Ha and \iPb) G 'Hb such that 

IV”) = IV'a) ® I'i/’s) ■ (1) 

In general, there are n := 2^“^ — 1 possible bipartitions of 'H and, if a decomposition of the form Q is 
found with respect to any one of them, then l"^) can be immediately declared biseparable. On the other 
hand, can only be declared GM-entangled after ruling out the existence of such a decomposition with 
respect to every bipartition of 'H. 

The notions of GM-entanglement and biseparability are extended to mixed states via a convex roof con¬ 
struction: An fV-partite mixed state p acting on 'H is GM-entangled if it is not biseparable; p is biseparable 
if it can be expressed as a convex combination of biseparable pure states, that is, 

P = ^Pi\'4’i){'<Pi\ ( 2 ) 

i 

where every is biseparable. Remarkably, although the biseparability of p requires every \'ipi) in Eq. 
to be biseparable, each one of them can be so with respect to a different bipartition of 'H, meaning that a 

biseparable mixed state need not to be separable with respect to any particular bipartition of 'H. 

So far we have only established the notion of GM-entanglement. In what follows we review some recent 
results on (i) how to tell if a given density matrix p is biseparable (detection) and (ii) if not, how to estimate 
the amount of GM-entanglement that it contains (quantification). Needless to say, these are hard problems 
even in the case N = 2 [33], let alone N > 2. 

2.2. GM-entanglement detection 

Detection of GM-entanglement has been intensely studied (see, e.g.. Ref. [38] for a review). To date, one 
of the most effective detection schemes was proposed by Huber et al. |33j . where it was shown that every 
biseparable state p satisfies every inequality of the |$)-parametrized family 

^|$)(p)<0, (3) 

where 

1$) - ^ ^(d>| (n^. ® IbJ^p^HUa, g IbJ |$) . (4) 

i=l 

In the above, |$) can be chosen as any product state of i.e., 

1^) = \Pi P 2 P 3 ■■■ Pn) G \ni 1^2 ^3 ■■■ vjs,). (5) 

The symbol II denotes the global permutation operator that peforms simultaneous permutations of all 
subsystems across the two copies of 'H, 

n(|/ri P 2 P 3 ■■■ Pn) ® \vi V 2 V 3 ■■■ vn)) = \vi V 2 V 3 ■■■ vn) ® \pi P 2 P 3 ■■■ Pn) , (6) 


whereas 11^^ G1^^ only permutes those subsystems whose labels are in Ai (of a given Hilbert space bipartition 
{Ai\Bi}). For example, consider the following bipartitions of the A^-partite Hilbert space 'H: 


Mil 

\Bi} = 

{1|2,3,. 


{A 2 \B, 

>.}-■ 

= {2|1,3,. 

..,N) 

and 

{Aat+i Hat+i} = {1,2 3, 


(7) 

Then, 














(Ha, 

G IbJ 

\PlP 2 P3 

■ ■ ■ Pn) 

G 

\viV2l'3 ■■ 

■ Vn) = 

= \vi P 2 

P3 ■ 

■ ■ Pn) G Ipi V 2 V 3 ■■ ■ 

Vn) , 

(8) 


(Ha, 

G ISa) 

\PlP 2 P3 

■ ■ ■ Pn) 

G 

WiV 2 V 3 ■■ 

■ Vn) = 

-- \Pl V 2 

P3 ■ 

■ ■ Pn) G \vi P 2 V 3 ■■■ 

Vn) , 

(9) 

(n 

Ajv+I G 

IfliV+l ) 

\PlP 2 P3 

■ ■ ■ Pn) 

G 

\vi V 2 V 3 ■■ 

• Vn) = 

-- \vi V 2 

P3 ■ ■ 

■ Pn) G \pi P 2 V 3 ■■ ■ 

Vn) ■ 

(10) 


Thus, if one can find a state |l>) such that (p) > 0, then p can be promptly declared to be GM-entangled. 
Unfortunately, such a criterion does not detect every GM-entanglement [there are GM-entangled density 
matrices that satisfy the entire family of inequalities ([^], but it is stronger than many commonly used 
criteria (see, e.g., m and references therein). 
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2.3. GM-concurrence estimation: The optimal-\^) scheme 

The relevance of transcends its application in the detection problem, manifesting itself also in 

the context of GM-entaglement quantification. Indeed, it was shown by Ma et al. m that, for every choice 
of 1$), the following inequality holds 


Cgme{p) > max [0,(p)] , 


( 11 ) 


where CgmeIp) denotes the concurrence of GM-entanglement of p; a well-defined (but difficult to compute) 
entanglement monotone. 

Unlike Cgme{p), the lower bounds 23\^){p) are easy to measure and to compute, emerging thus as 
natural candidates for GM-concurrence estimation. However, the accuracy of such estimates are strongly 
dependent on the choice of |$) for a given p, highlighting the importance of devising efficient ways of making 
such a choice. 


Formally, the strongest inequality of the family (111 is given by 


Cgme{p) > max 


0,2rmix.?l$)(p) 


= C'|$)(P): 


( 12 ) 


where the optimization runs over all product states of Since any product state can be constructed 

from LU transformations on a reference product state |$o) G we can set 


( N N \ 

(g)KG(g)W^n po) , 

n=l n=l / 


(13) 


where Vn and Wn are SU{dn) elements. In this case, the optimization is performed over 2N independent 
(special) unitary transformations (Vn and Wn), in terms of which the objective function takes the form 


^l$)(p) = I (0| ulpUo |0) I - ^7(0|«7/pi7,|0) {0\ulpU,\0), 


(14) 


2=1 


where we have chosen |<i)o) = |0)®^ G |0)®'^ = |0) G |0) and defined the LU transformations 


N 

N 

N 

N 


Gs 

II 

Uo := (g) , 

n=l 

Ui := (^) Si^n , 
n=l 

U, :=(g)S,,„. 

n—1 

(15) 


In the above, the single particle unitary matrices Si^n and Si^n have the index i referring to a particular 
bipartition {Ai\Bi} of TL and, in terms of this, are given by: 


Wn if n e 
Vn otherwise 


and Si^n 


Vn if n e 
Wn otherwise 


(16) 


Thus, we end up with the optimization problem 


max 

v„,w^esu{d„) 


{0\ulpUo\0) 


j2\/{0\UlpU,\0) {0\ulpU,\0). 

2=1 


(17) 


Let us now restrict the analysis of this optimization problem to the case of iV-qubits (i.e., dn = 2 for 
every n = 1,...,N). Since single-qubit unitary transformation require only 2 real parameters and our 
objective function involves 2N unitaries, the number of real variables to be optimized is 4iV. Although the 
linear scaling with N is certainly an appealing feature of this approach, the resulting objective function is 
not everywhere differentiable. This fact compromises the performance of numerical algorithms that rely on 
differentiating the objective function to approach the optima. 
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3. X-Heuristic 


We now introduce a heuristic approach for estimating the GM-concurrence of a density matrix of N- 
qubits. To motivate our heuristic, we start by reviewing the concept of X-density matrices and how to 
quantify their GM-concurrence. An A^-qubit A-density matrix is any density matrix that, in the computa¬ 
tional basis, takes the form 


Px = 


ai 


02 


^26 


i<i>2 


rie 


iipi 






r2e 




(18) 


with n = 2-^“^ and, for every fc = 1,... ,n, {ak,hk} G R-i-, (t>k G [0, 27r], G [0, y/akbk] (non-negativity of 
Px) and + ^fc) = 1 (normalization of px)- As shown by Hashemi Rafsanjani et al. [22], such states 

can have their GM-concurrence easily — and exactly — quantified by the formula: 


Cxipx) = max[0, cx(px)] 


(19) 


where 


cx{px) '■= 2 max 

feGe„ 


Tk 


- 


jee„\{k} 


and &n ■= {1, ■ ■ •, ■ 


( 20 ) 


The proposed heuristic takes advantage of the simplicity of Eq. (19) to produce computable estimates of 
GM-concurrence for non-X density matrices. In fact, as noted in Ref. [22] (and explictly demonstrated 
in the Appendix A), the evaluation of Cx{p) for an arbitrary A-qubit density matrix p (written in the 


computational basis), leads to a lower bound on the GM-concurrence of p. In principle, depending on the 
particular form of p, this lower bound can be very loose and, consequently, a very poor estimate of Cgme{p) 
[see, e.g., Eqs. ( A.15| ) and (A.16)]. However, since this lower bound is saturated for X-density matrices, one 
should expect that the “closest” p is to a X-density matrix, the better that estimate will turn out to be. 

The X-heuristic materializes this reasoning by “approximating” p to a X-matrix with a LU transformation 
Ui and, subsequently, estimating its GM-concurrence as Cx{p), where 


p.= UipUi 


t 


( 21 ) 


Of course, the constraint of LU transformation is sufficient to imply that p and p have the same amount of 
GM-concurrence. However, whereas Cgme{p) is invariant under LU transformations of p, Cx{p) is not. 

Once again, we are left with the problem of choosing a LU transformation, but now aiming to implement 
the “approximation to the X-form”. In order to formalize the meaning of this approximation, we introduce 
the penalty function 

fip)= E (22) 

j>i 

j+z#2~ + l 

which (i) vanishes if and only if p is a X-density matrix and (ii) smoothly increases as the magnitudes of 
the off-X entries of p increase. Under this “/-metric”]^ the optimal choice of Ui is given by the solution of 
the optimization problem 

rmnfiUipU^), (23) 

Ui 


^Here and throughout, the term “metric” is used in a loose way to suggest a notion of distance between a given density 
matrix and the set of X-density matrices. 
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which is a non-linear least square problem on 2N real variables with residuals given by \pi,j\- We can think 
of it as finding the LU related density matrix that most accurately “fits” the X-density matrix model. 

Besides halving the number of variables in the optimization problem 0 . the objective function in ( [ 2 ^ 
is smooth, which implies that numerical algorithms can take advantage of well-defined derivatives to more 
efficiently approach a minimum. On the analytical side, for density matrices p with enough symmetry, the 
objective function fCUipUj) will be sufficiently simple to enable educated guesses on the location of the 
minima, which can be certified a posteriori by evaluating the gradient and the spectrum of the Hessian 
matrix at the candidate point. In any case, even if the gradient vanishes and the Hessian turns out to 
be positive definite, we cannot be sure that we have attained a global minimum. Nevertheless, at this 
point we may have already obtained a non-trivial lower bound on the amount of GM-concurrence of p at a 
considerably low computational cost. 


4. Applications 

In this section we apply the X-heuristic to estimate the amount of GM-concurrence in two families of 
A-qubit density matrices: the diagonal symmetric states and the steady states of the driven Dicke Model 
at zero temperature. For each family, we benchmark the quality of the resulting estimates against those 
produced by the optimal- 1 $) scheme. 


4.I. Diagonal symmetric states 

Diagonal symmetric states are a natural extension of W-states to the domain of multiqubit mixed states 
and have recently attracted much attention in the helds of GM-entanglement detection and quantification [331 
1331 [S3HS3] • Formally, they can be defined as 

N 

Pds,N:=J2Pk\Dk){D^\ , ( 24 ) 

k=0 


where pk G [0,1] and J2kPk = I- Here, the states \D^) are the (totally symmetric) A-qubit Dicke states of 
k excitations (321 [SS] , defined in the computational basis as 


D^) := 


0 k(i, 


-•,l,0,^-^0)) 


(25) 


with the summation running over every distinct permutation of the sequence of k ones and N — k zeros. Of 
course, there are Cj^ = A!/[fc!(A — fc)!] such permutations, justifying the normalization constant upfront. 
Moreover, note that \D^) are the A-qubit W-states |S5], which sets the context in which pds,N can be 
considered a generalization thereof. In the Appendix B we give explicit matrix forms of Pds,N in tbe 
computational basis for A = 2, 3,4. 

For A = 2, the resulting density matrix is already an X-state [cf. Eq. (B.l)]. As a result, Ui = I 4 , 
fmin — 0 and CxiPdsp) = CoMEiPdsp)- For A > 2, the corresponding density matrices are no longer of 
the X-form [see, e.g., Eqs. (B.2| and (B.3)] and thus, in order to optimally estimate Come with Cx, we 
must first determine a LU transformation 


Ui 



cos'dj sini?je*'^J 

— siniJje”*'^^ cos'dj 


(26) 


that minimizes f{UiPdB,NUj). 
reveals a minimum ai[ 3 ] 


Thanks to the symmetry of pds,N^ some simple analysis of this function 

TT 

dj = — and (pj = 0 for every j € ©at , (27) 


complementary numerical analysis strongly suggests that Eq. (|27|l gives, actually, 


a global minimum. 


6 











for which the corresponding values of / = /mm can be promptly computed with the aid of Eqs. (21) and 
(221. In Fig. these are plotted for a family of populations pk specified by the single parameter r € [0,1] 
as follows: 



Figure 1: Distance (“/-metric”) between the set of X-density matrices and the LU transformed diagonal symmetric state as a 
function of r. The included dashed (red) lines represent (half of) the sum of the absolute values squared of the anti-diagonal 
terms. For N = 2, the density matrix is of the X-form for every value of r E [0,1] and, as r increases, it approximates a 
diagonal matrix. For N > 2 the density matrices are not of the X-form for any value of r E [0,1] and better approximate this 
condition at intermediate values of r. 


Intuitively, the plots of Fig. can be interpreted as some error measure associated with the process of 
estimating CGME{Pds,N) as Cx{Pds,N)- Indeed, if /mm = 0 (e.g., case N = 2), then Cx{pds,N) is exactly 
equal to CQME{Pds,N)- As / increases, pds,N deviates from the X-form and, hence, the estimate CxiPds,N) 
is expected to be poorer. Nevertheless, recall that Cx(pds,N) is always a lower bound on CGMEiPds,N), so 
the “/-error bars” are not centered in Cx{pds,N), but lay strictly above it. 

In Fig. we present the resulting values of Cx{Pds,N) (solid line) along with the points corresponding 
to Ci^) (pds.Tv) for one hundred values of r uniformly distributed in the range [0,1]. While the line was 
analytically constructed, each point was obtained by numerically solving the optimization problem (Hill- 
Overall, we notice that both estimates follow very similar trends and coincide in a significant portion of the 
domain. Next, their similarities and differences are discussed in greater detail. 

For N = 2, owing to the X-form of Pds, 2 , Cx is precisely equal to the Wootters concurrence and, therefore, 
exactly matches for every r G [0,1] — see Plotj^. For N = 4,6, this exact matching no longer occurs 
for 0 < r < 0.212 {N = 4) and 0 < r < 0.I2I {N = 6). Nevertheless, as the insets in Plots [^,c show, the 
difference between and Cx is orders of magnitude smaller than the actual value of the estimates, being 
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a N = 2 


b TV = 4 


c iV = 6 



0 1 2 3 4 1 

5 5 5 5 




e N = h f N = 1 




Figure 2: Comparison between lower bounds (LBs) on the GM-concurrence of a family of A^-qubit diagonal symmetric states 
[cf. Eq. (|28|l]. The solid line represents Cx{pds,N) the red circles represent Cj^{pds,N)- While the density matrices 
Pds,N (and their corresponding X-concurrence) were analitically computed [cf. Eqs. | |24[ |-( [2^ ]) the values of C\^^^{pds,N) were 
obtained from a numerical implementation of the optimal |$} method [cf. Eq. 


thus negligible for most practical applications. For odd V, a more significant mismatch between Cx and 
occurs at the extremes of the domain (cf. Plots [^,e,f), but it rapidly decreases and vanishes in a large 
central region of the domain, where neither estimates detect GM-entanglement. An important difference 
between Cx{pds,N) and {pds,N) becomes apparent in the case N = 7 (Plotj^), where we note that each 
estimate arrives to (and departs from) zero at slightly different values of r, in such a way that there are a 
few states whose GM-entanglement can be detected by the optimal-1$) scheme, but not by the X-heuristic. 

4-2. Steady states of the driven Dicke model 

From the theoretical point of view, the Dicke model |43j represents a special quantum mechanical model 
whose unique mathematical and physical properties allow us to describe important cooperative phenomena 
(such as, for example, resonance fluorescence and cooperative emission) far from the thermodynamic limit. 
Their virtues are also reflected in works involving the exact steady state solutions of its master equation 
that include (or not) the effects of detuning between a collective driving field and the atomic resonant 
frequency [46l-j49]. In this section, we focus on the steady state solution of the driven Dicke model at zero 
temperature 


Ps,N 



(29) 


where N depicts the number of two-level atoms (ions), Dx is a normalization constant, 7 := 7 ^/D denotes 
the ratio between the Einstein A-coefficient 7 ^ of each atom (ion) and the Rabi frequency D, and g '■= if j. 
Note that J± := G crj- G l 2 w-i for ^ = 1 ,..., V correspond to the collective raising and lowering 
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operators expressed in terms of the Pauli matrices (t± which satisfy, together with the collective inversion 
operator the usual angular momentum commutation relations. 

In the Appendix B[ explicit matrix forms of Ps,N in the computational basis are given for TV = 2,3,4. 
Noticeably, the resulting matrices are not of the X-form and, as shown in Fig. they cannot be exactly 
brought to the X-form via LU transformations (i.e., they do not assume the X-form in any orthonormal 
basis of product states). These observations qualify ps jv as a good testbed for the X-heuristic and so, in 
Fig. HI we present the corresponding estimates obtained from numerical implementations of the X-heuristic 
and optimal -1 $) scheme. 









0 1 


•10 


-3 


2 3 4 5 6 7 

= 5 


0 1 


5 6 7 



•10 


-3 


= 6 



•10 


-3 


N = 7 



Figure 3: Distance (“/-metric”), as function of 7 , between the set of X-density matrices and the closest LU transformed steady 
state of the driven Dicke model. The included red data points represent (half of) the sum of the absolute values squared of 
the anti-diagonal terms. For N > 2, the LU operation that optimally decreases the off-X entries has a deleterious effect on the 
anti-diagonal as well, essentially making diagonal the transformed density matrix. 


At a glance of Fig. both Cx and are seen to follow very similar trends: a zero plateau that extends 
from 7 = 0 up until 7 « 1, followed by a sudden growth and, finally, an asymptotic decay. In what follows 
we offer a more in-depth analysis of these results by considering the cases N = 2 and N > 2 separately. 

In the case A^ = 2 , on top of Cx and , we have also plotted the Wootters concurrence of Ps ,2 which, 
in this case, is the “genuine multipartite” concurrence Cgme- Despite ps,2 not being of the X-form, we 
note that the three measures coincide up to the numerical precision of 10 “^'^, with the zero plateu ranging 


up until 7 = 1.00 and the maximal concurrence (~ 7.735 x 10 occurring at 7 „ 


1.65. Interestingly, 


the observed coincidence between the three measures signals that the X-concurrence formula may also be 
an exact concurrence formula for some non-X states. In [Appendix A this point is further explored with 
a demonstration that Cx{Ps, 2 ) matches exactly the Wootters concurrence (Cw) of Ps^ 2 - This observation 
provides a (constructive) negative answer to an open question in Ref. |50j on whether the saturation of 
Cw{p) > C'x(p) requires p to be of the X-form in some orthonormal product basis. 

For N > 2, Cx and no longer coincide for every value of 7 . Nevertheless, the difference between 
them asymptoticaly decreases and nearly vanishes at relatively small values of 7 . For example, the ratio 
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Figure 4: Comparison between the lower bounds Cx{Ps,n) (black circles) and ( Ps.n ) ('"6^ circles) of the genuine multi¬ 
partite concurrence of the A'^-qubit steady states of the driven Dicke model [cf. Eq. l |2’9[ |]. Both estimates were numerically 
calculated. 


— Cx)/C\^'f is smaller than 3.2% for 7 > 7 (iV = 3), 7 > 2.2 {N = 4), 7 > 1.66 {N = 5), 7 > 1.4 
(TV = 6) and 7 > 1.2 (TV = 7). Of course, this behavior was to be expected since for such values of 7 we 
have already seen, in Fig. that Ps,N is nearly X-formed. In fact, thanks to that, we can also infer that 
Cx must well aproximate the (unknown) value of Come in such values of 7 . 

As it ocurred with the diagonal symmetric states, the main drawback in replacing with Cx is that 
there are certain states that have its GM-entangled detected by the former but not by the latter. In Fig. 
this is expressed by the fact that the zero plateau generated by the X-heuristic is longer than that obtained 
with the optimal-1 $) scheme. However, as the plots show, the difference is usually very small, meaning 
that failure to detect GM-entanglement with the X-heuristic will only occur for very few and specific values 
of 7 . Other than that, it is also worth noticing that each method foresees a maximal value of Come at 
slightly different values of 7 , in such a way that 7 ^L > 7 max . However, the difference — 7 max is 

again very small and, ultimately, it is not clear which one more accurately describes the actual location of 
the maximum of Cqme- 

5. Computational Efficiency 

The GM-concurrence estimates produced by the X-heuristic cannot improve on the corresponding esti¬ 
mates produced by the optimal-14)) scheme; a fact that follows directly from the derivation of the X-heuristic 
(cf. Sec.[^. There is, however, an important trade-off between accuracy and efficiency that must be consid¬ 
ered before disregarding the X-heuristic in favor of the optimal-14>) scheme. On the accuracy side, we have 
already seen that the X-estimates are usually pretty good approximations of the optimal-|4>) estimates. On 
the efficiency side, we now explicitly demonstrate that the X-estimates are significantly easier to compute 
both for the diagonal symmetric states and for the steady states of the driven Dicke model. 
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As mentioned before, the computational advantage of the X-heuristic over the optimal-1$) scheme relies 
upon two aspects of the objective function associated with the former: (i) it is smooth and (ii) it depends 
on a number of variables that halves the number of variables in the (non-smooth) objective function of the 
optimal-1$) scheme. While the relevance of (ii) is obvious, the importance of (i) cannot be overestimated. 
Thanks to (i), the X-heuristic optimization problem can be regarded as ideal for the application of a quasi- 
Newton method a family of nonlinear optimization algorithms that take advantage of well-defined 
derivatives of smooth objective functions to more efficiently converge to a minimum. 

In this work, the numerical computations of the X-heuristic and optimal-|<i>) estimates were performed 
with the MATLAB function fminunc, which implements a popular quasi-Newton algorithm known as BFGS 
(after its discoverers Broyden, Fletcher, Goldfarb and Shanno). As with any quasi-Newton algorithm, the 
BFGS performs successive evaluations of the gradient to build a quadratic model of the objective function 
that is sufhciently good to attain a superlinear rate of convergence. It contrasts with Newton’s method m in 
the sense that it does not require (nor attempts to compute) the Hessian matrix, typically a time-consuming 
and error-prone task. Instead, at each step, it gains information about the second derivative along the 
search direction by considering changes in the gradient. Although the BFGS algorithm requires repeated 
computations of the gradient, it has been noted to perform well also in nonsmooth optimization problems, as 
long as it does not run into a nonsmooth point (see, e.g., [5S] and references therein). Consistently with this 
observation, we have observed a better performance of the optimal-1$) scheme with the BFGS algorithm 
than with algorithms that do not rely on evaluations of the gradient (e.g., Nelder Mead algorithm [59] . 
implemented by the MATLAB function fminsearch), especially for larger values of N. For this reason, 
aiming to build the most efficient implementation of the two schemes, we employed the BFGS algorithm for 
both. 

In order to compare the performance of the two schemes in producing reasonable estimates of GM- 
concurrence, we set a threshold equal to the X-estimates presented in Figs. and and timed how long it 
took for each scheme to reach that threshold. As it is usually the case with search algorithms, a starting 
value had to be provided, and so we resorted to random initial guesses (uniformly sampled) to avoid biasing 
the optimizer toward (or against) a satisfactory minimum. By doing so, we could gauge how hard it is for 
each scheme to reach the threshold from a zero-knowledge initial condition. 

Each optimization was repeated 100 time^ and the five-number summary of the resulting “wall-clock 
time” distributions is shown in Fig. [^as a box-and-whiskers plot. Plot[^ refers to the time measurements 
obtained for the diagonal symmetric state with r = 3.03 x 10“^ for every N G [2, 7], whereas Plot IJr refers 
to the time measurements for the steady states of the driven Dicke model that, for each N, have the value 
of 7 corresponding to the peak of the X-estimates in Fig. We note that to obtain each box-and-whisker 
appearing in Fig. several initial guesses were typically made because either (i) the optimizer converged 
to a unsatisfactory minimum or (ii) the optimizer failed to converge after 10"^ iterations. Of course, the 
resulting box-and-whiskers take into account every (if any) unsuccessful attempt as well as the successful 
one. 

Remarkably, Fig. reveals a significant advantage of the X-heuristic over the optimal-1$) scheme, es¬ 
pecially for the larger values of N. For N = 7, for example, the X-heuristic estimates were produced in 
roughly 10 s, whereas the optimal-1$) scheme was unable to reach the threshold once, even after a week 
(approximately 6 x 10® s) of computation. Most importantly, the plots suggest that the computational time 
for each scheme scales differently with N, with the X-heuristic resembling an exponential growth and the 
optimal-1$) some worse scaling. Naturally, this observation cannot be explained by the difference in the 
number of variables associated to each scheme (since they both scale linearly with N), and we assign it to 
the smoothness (or lack thereof) of the objective functions. Because of this scaling, the X-heuristic may be 
the only viable option for GM-concurrence estimation in generic A^-qubit states with intermediate values of 
N. 

The results of this section indicate that, despite the theoretical superiority of the optimal-|$) scheme 
over the X-heuristic, the latter may actually overperform the former in practical situations where a time 


^Except for the time measurements of the optimal-|<I>) scheme for states with N = 6, which were repeated only 5 times due 
to the long time-frame necessary to converge to a satisfactory value. 
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Figure 5: Box-and-whisker plots of computation time taken to reach the X-estimates of GM-concurrence shown in Figs. and 
l^for (a) the A^-qubit diagonal symmetric states with r = 3.03 X 10“^ and (b) the A'^-qubit steady state of the driven Dicke 
model (zero temperature) corresponding to the peaks of the X-estimate in Fig.|^ i.e., 7 = 1.652 (N = 2), 7 = 1.623 {N = 3), 
7 = 1.362 {N = 4), 7 = 1.217 {N = 5), 7 = 1.072 {N = 6 ), 7 = 0.956 {N = 7). The time measurements were conducted while 
running our MATLAB implementations of the X-heuristic (black) and optimal-|'l>} scheme (red) with random initial guesses 
(uniformly sampled) and termination tolerances on the parameter values (TolX) and on the objective function value (TolFun) 
set to 10“^^. For each combination of method and state, the timing was repeated 100 times (except for the optimal-|4>) method 
applied to states with A = 6 , in which cases we contented ourselves with only 5 repetitions). The plots display the five-number 
summary of the resulting time distributions. Whenever the optimization terminated without converging to a value equal to 
or greater than the desired threshold, another try was made (with a different random initial guess) until the threshold was 
reached. Computations were performed on a 2.2-GHz Intel Core i7-2670QM. 


horizon has to be considered. For exampie, if we aiiow one second to obtain a GM-concurrence estimate for 
a state with N G [3, 5], then Fig. [^impiies that we are much more iikeiy to obtain a better estimate with the 
X-heuristic than with the optimai-|$) scheme, which may have converged to a iesser estimate (or even have 
not converged at aii). Naturaiiy, time constraints iike this become virtuaiiy more important as TV increases. 

6. Concluding Remarks 

In this paper we introduced a heuristic method for estimating the GM-concurrence of an A^-qubit density 
matrix. It consists of evaluating the TV-qubit X-state GM-concurrence formula for density matrices that are 
not of the X-form, but were previously brought “as close as possible” to it with a LU transformation. We 
have shown that the estimates thus produced are lower bounds on the GM-concurrence by demonstrating 
that they actually lower bound a standard lower bound on the GM-concurrence (attainable via a numerical 
procedure which we refer to as the “optimal-1$)” method). Most importantly, by examining two prominent 
families of mixed A^-qubit states, we found that our estimates are usually very close to those produced by 
the optimal-1$) scheme and are significantly easier to be numerically computed and, in certain cases of high 
symmetry, can be analytically obtained. 

As for future directions, let us mention two possible extensions of the present work. First, it would be 
interesting to evaluate the performance of the X-heuristic in other well-established families of A^-qubit mixed 
states such as reduced ground states of spin chain models. As a matter of fact, it is fair to say that the 
X-heuristic has already been successfully applied in this context, in Ref. [21], where the genuine tripartite 
concurrence of the reduced ground state of three spins symmetrically distributed in the cluster-Ising model 
was exactly calculated thanks to the fact that such states can be turned into X-states via LU transformations. 
Nevertheless, the application of our heuristic in alike models for which the X-form is not exactly attainable 
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is yet to be explored and is likely to provide an efficient tool for GM-entanglement estimation in condensed 
matter systems. 

Second, the requirement of approximating the X-form via LU transformations is unnecessarily strong for 
the goal of preserving GM-concurrence. In fact, it was recently conjectured that a generic two-qubit density 
matrix can always be turned into a X-density matrix of same entanglement via a unitary transformation 
that is not necessarily local [SD] • In Ref. EH, this conjecture was rigourously proved for three entanglement 
measures and, moreover, a semi-analitic prescription for constructing the corresponding “entanglement- 
preserving-unitaries” was delineated. Along the same lines, it would be interesting to characterize the 
most general family of “GM-concurrence-preserving-unitaries” and attempt to approximate the X-form via 
conjugation with elements of this less constrained family. In that case, better X-form approximations should 
be expected, presumably improving the quality of the resulting GM-concurrence estimates. Of course, this 
new optimization problem may be less suitable for numerical implementation, in which case the trade-off 
between efficiency and accuracy should be carefully reconsidered. 


Acknowledgements 

P.E.M.F.M. acknowledges the financial support of the Brazilian Air Force and the program “Ciencia sem 
Fronteiras”, Project No. 200024/2014-0. 


Appendix A. X-concurrence is a lower bound for GM-concurrence 


In this appendix we demonstrate that the X-concurrence formula (19) applied to a generic A^-qubit 


density matrix p is equal to max^ (p) ^ certain family of states {|4 ’^)}/Zq (recall that n := 2"“^), 

being thus a lower bound on the GM-concurrence. 

Consider the following (decimal) representation for the A-qubit computational basis 


|0) = |00 ••• 00) , |1) = |00 ••• 01) , ... , |A) = |10 ••• 00) , 

|A + 1) = |10---01) ,..., |2'^-1) = |11 ••• 11) , (A.l) 

and, in terms of this, define 

1$^) := l/r) G 12-^ — 1 — ^) for every /r e {0,..., n — 1} . (A.2) 

Substitution of this into Fq. gives, after some straightforward algebra, 

= \J {'^\ P k) (2^ - 1 - i^l p |2^ - 1 - u). (A.3) 


Now, let the density matrix of p be decomposed (in the computational basis) as 


P = 


Oi 


02 


r 2 e 


*02 


rie*' 


On ^nC 

g-i0„ 5 


i(l)n 


r2e 




Tie 


o 


(A.4) 


where O is some 2^ x 2^ Hermitian matrix with zeros along the main- and anti-diagonals. In terms of this, 
Fq. (A.3) clearly takes the form 


■^\'S>p,){p) —- 'Y a / a^+iK+i- 

!/G{0,...,n-l}\{/i} 


(A.5) 
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Making i/ + 1 = j and ^ + 1 = fc (in such a way that j, k G S„), it is clear that 


and, thus [cf. Eq. (20)] 


= n- \/ajVj 

jee^\{k} 


cx{p) = 2 max . 

kG&n 


(A.6) 

(A.7) 


This fact, combined with inequality (12) and identity (19), leads to the desired result: 
Cgme[p) > max 


0 , 2 nuixJ?l$)(p) 


> max 


0,2 rnax 

kG&n 


= max[ 0 , cx (p)] = Cx{p). (A. 8 ) 


As suggested in Ref. [S^ (within the framework of TV = 2 ), it is interesting to look for density matrices 
that saturate the above inequalities. Of course, saturation occurs if p is an X-state, other than that (and to 
the best of our knowledge), no examples of density matrices saturating inequality (A. 8 ) have been identified 
so far. In what follows, we introduce a family of two-qubit density matrices which are not of the X-form 
in any product state basis and, nonetheless, satisfy Cw{p) = Cx{p), where Cw denotes the Wootters 
concurrence (i.e., the “genuine multipartite” concurrence of two-qubit states). 

Consider the two-qubit steady states of the driven Dicke model at zero temperature, whose density 
matrix in the computational basis is presented in Eq. (B. 8 ). Noticeably, it is not an X-state for 7 7 ^ 0 and, 
as Plot ii shows, there is no LU transformation capable of simulatenously vanishing all of its non-X entries 
(although this can be remarkably well aproximated for sufficiently large values of 7 ). Nevertheless, as the 
following computation shows, the Wootters concurrence exactly matches the X-concurrence in this case. 

We first apply the X-concurrence formula to Eq. (B. 8 ), which yields 


Cx{Ps, 2 ) = 2 max 


- 1 7^ - -I- 272 -I- 474 


Do 


Do 


= max 


2 ( 7 ^- 1 ) 

Do 


(A.9) 


Then, in order to compute the Wootters concurrence of Ps, 2 : we start by evaluating 


P^. 2 (« 


>)[ps,2nc 


v ) — 7 -i 2 


Dl 


2^73 2173 -272 

4173 1 - 274 -274 -2I73 

4^73 - 27 ^ 1 - 27 ^ -2173 

-272 - 87® -4173 -4173 1 -H 47^ 


(A.IO) 


whose eigenvalues are 


1-^274-^272^1-^74 1 1 1-^274-272^^1-^74' 


D 2 

^2 


’Dr Dr 


n 2 

^2 


Since the first element in the above set is clearly the largest eigenvalue, one has 


Cw{Ps, 2 ) = max 

^2 


= max 


0, 71-^274-^272^1-^74 - 2 - 7 i + 27 ^- 27 V 1 + 7 ^ 

2(7^ - 1 ) 

Do 


where, in order to obtain the second row, we used the denesting identities 

7l + 274 ± 272 a/T-I- 74 = vTTqd ± 72 ^ 


(A.ll) 

(A.12) 

(A.13) 

(A.14) 


Comparing Eqs. (A.9) and (A.13), it is clear that Cw{ps, 2 ) = Cx{Ps, 2 )- 
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Unfortunately, saturation of (A. 8 ) does not occur for the zero-temperature steady states of the driven 
Dicke model with N > 2. For N = 3,4, for example, density matrices in the computational basis are 
explicitly presented in Eqs. (B.9) and (B.IO), respectively. Although Plots |^,c reveal that these states are 
GM-entangled for 7 > 1 , a straightforward computation of theirs X-concurrence yields 


Cx{Ps,3) = max [ 0 , 673 - 3yri(r2 +474), 27 ^ - 2v/r\(f^+^ - + 

L>3 L 


= 0 (A.15) 


and 


Cx{Psa) = yw max [o, 247 ^ - 3 (r 2 + 47 ^) - + , 

IJ4 L 

67 ^-3(r2 + 47^) -3x/rdryrTT^- \/r7T24^W+5^, 

47^ - 2(r2 + 47^) - 4 a/UJVTTT27^ - = 0 , (A.16) 


which shows that Cx is only a trivial lower bound for Cqme in these cases. 


Appendix B. Density Matrices in the Computational Basis 

In this appendix we present the density matrices, written in the computational basis, of the A-qubit 
diagonal symmetric states and steady states of the driven Dicke model, for N = 2,3,4. For ease of visualuza- 
tion, we restrict to show the lower triangular entries, replace zeros with dots and boldify the anti-diagonal 
entries. 

Appendix B.l. Diagonal Symmetric States 
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Pds,2 — 2 
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(B.l) 
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(B.2) 


15 






















Pds,A = 


4po 



Pi 












Pi 

Pi 

P2 










Pi 

Pi 

P2 

Pi 

P2 










P2 


P2 

P2 






1 







P3 





4 

Pi 

Pi 

P2 

Pi 

P2 

P2 

• Pi 

P2 







P2 

. 

P2 

P2 

P3 ■ 

P2 

P2 

Pz 





P2 


P2 

P2 


P2 

P2 


P2 



• 





Ps ■ 



Pz 

■ Pz 


• 






Pz ■ 



Pz 

■ Pz 


4p4 


(B.3) 


Appendix B.2. Steady states of the driven Dicke model 

For brevity, we define the following family (parametrized by k) of quadratic functions of 7 : 


Ffc := 1 + k-f^ . (B. 4 ) 

The normalization appearing as a multiplicative factor in front of the matrices forms are given below: 

:= 4(1 + 72 + 7^), (B. 5 ) 

D3 := 4(2 + 37 ^ + 67'* + 97®), (B.6) 

£>4 := 16(1 + 272 + 67^ + I87® + 367®). (B. 7 ) 
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1 + 7^ 
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7^ 

1+7^ 


-272 

17(1 + 27^) 

17(1 + 27^) 1 + 27^+47^ 
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'yir2 7^r4 7^r4 —227^r6 47^+r2 
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